Spin-flip induced magnetoresistance in positionally disordered organic solids 
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A model for magnetoresistance in positionally disordered organic materials is presented and solved using 
percolation theory. The model describes the effects of spin flips on hopping transport by considering the effect 
of spin dynamics on an effective density of hopping sites. Faster spin-flip transitions open up 'spin-blocked' 
pathways to become viable conduction channels and hence produces magnetoresistance. The magnetoresistance 
can be found analytically in several regimes, including when the spin-flip time is slower than the hopping time. 
The ratio of hopping time to the hyperfine precession time is a crucial quantity in determining the shape of 
magnetoresistance curves. Studies of magnetoresistance in known systems with controllable positional disorder 
would provide a stringent test of this model. 



Spintronics 1 1 ] in organic materials has generated consider- 
able interest in recent years |2], due to the long spin lifetimes 
of organic semiconductors as well as the flexibility, low cost 
and chemical tunability of organic devices |3|. Spin transport 
properties are intimately connected to the electrical transport 
properties |4], so although spin transport through inorganic 
semiconductors has been extensively explored O, novel 
features should be expected in organics due to their very dif- 
ferent electronic transport properties. The understanding of 
spin transport in organics has been challenged by the discov- 
ery of magnetic field effects on properties such as conductivity 
and electroluminescence I7UT511; and characterized by magne- 
toresistances of 10-20% in magnetic fields as small as 10 mT. 
Several new models of organic magnetoresistance (OMAR) 
have been proposed, many of which involve spin-dependent 
processes emanating from hyperfine interactions. These mod- 
els can be broadly categorized by their reliance on bipolaron 
lfl~3ll or electron-hole pair J51H01 formation. Their main points 
of difference are that bipolaron models consider the relative 
spins of like charge carriers and electron-hole models con- 
sider the relative spins of electrons and holes on different sites 
in the formation probability for a two-carrier entity. However, 
no model of OMAR has explicitly taken into account how the 
presence of spin-blocked sites affects the theoretical descrip- 
tion of hopping transport for a single carrier using percolation 
theory lfT6lfT7l 

This Letter provides a description of magnetoresistance 
based on percolative hopping transport for positionally dis- 
ordered organic semiconductors. The model proposed here 
maps the complex phenomena of spin-dependent hopping 
onto a simple problem of r-percolation with an effective den- 
sity of hopping-accessible sites that depends on the magnetic 
field through spin relaxation. We focus on unipolar charge 
transport since several analytic results can be readily obtained; 
extension to bipolar transport can be done with similar tech- 
niques. Our percolation based model allows us to explain the 
width and saturation of measured magnetoresistance curves 
as well as make predictions of magnetoresistance in systems 
with low site concentration and high temperatures. Finally, 
we propose experiments that could test our theory and thereby 
shed light on the operative mechanisms leading to OMAR. 

Model - We model the spatially disordered organic system 



as a network of random resistors in the spirit of Miller and 
Abrahams 1 18]. The resistance between two sites, i and j, is 
given by Rij = R^e lri ^ a where is the their separation and a 
is the localization length of a carrier at a site which we assume 
to be constant throughout the system. The bulk resistance in 
such a random resistor network is solved by percolation the- 
ory (l6l[T9][20]|. The bulk resistance is governed by a critical 
resistance (distance) R c (r c ) which is the smallest resistance 
(or equivalently the smallest separation) that still allows for 
an infinitely large network of bonds. This percolation length 
is set by the bonding criterion: B c = 4% r 2 Ndr, where N 
is the density of sites in the system and B c is a number that 
determines how many bonds each site in the percolating net- 
work must connect to on average; B c « 2.7 in three dimen- 
sions |[T6l . Energy disorder is negligible when the inter-site 
separation is large and temperatures are high; conductivity 
due to r-percolation has been observed in organic semicon- 
ductors in this regime (SUE). These conditions are assumed 
throughout this Letter. For smaller inter-site separations or 
lower temperatures, energy disorder plays a pivotal role. Our 
model could, in principle, be generalized and solved numeri- 
cally to treat such situations. 

Spin affects electronic transport in hopping transport 
through the Pauli exclusion principle, as shown schematically 
in Fig[T} double occupation of carriers on a single site is for- 
bidden if their spins are parallel (P), but allowed if they are 
anti-parallel (AP) [Q21 E3 (at the cost of a Coulomb interac- 
tion energy U). To clarify the role of spin blocking we use 
the simple picture of U = 0, so a carrier with arbitrary spin 
is restricted from hopping to an occupied site with P-spin but 
may hop to a site occupied by an AP-spin just as it would 
to an unoccupied site. The respective concentrations of these 
three types of sites are Np, Nap, and No. We consider carrier 
concentrations dilute enough to neglect hops to doubly occu- 
pied sites. Since carrier hopping to an occupied site with a 
P-spin is forbidden, the concentration of sites is effectively 
reduced to N — Np. In the absence of spin flips we would 
then write the bonding criterion as B c = 4n J^ c r 2 N f e ^dr where 
N' e ff = N — Np. The spin flip of a carrier at a site can be un- 
derstood as a dynamical process that will cause the relative 
spin orientation between two singly occupied sites to change. 
Hence, the hopping dynamics between two occupied sites is 
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FIG. 1. Spin blocking in transport. Top: the initial situation for 
hopping when a carrier's nearest neighbor is occupied by a parallel 
spin (left) and anti-parallel spin (right). Bottom: on the left the spin- 
blocked carrier has made the more difficult further hop to an unoccu- 
pied site. On the right, the carrier successfully hops to the occupied 
nearest neighbor with anti-parallel spin. The carrier concentration is 
dilute enough such that if a carrier's nearest neighbor is occupied, it 
is improbable that its second nearest neighbor will be occupied. 



strongly dependent on spin flips. If the total concentration of 
singly occupied sites is fixed at N s , then at any given time the 
average densities of P-spin and AP-spin sites are N s /2. Thus 
as a carrier attempts a hop to a singly occupied site, the prob- 
ability of success will be 1/2, independent of spin effects. So, 
half the time the hop will be successful and the density of sites 
for which these successful hops take place is Nap = N s /2. So 
as before the density of unrestricted hopping sites is N' e ^. We 
must now account for the situation that occurs the other one- 
half time in which the hopping attempt is foiled due to the 
singly occupied site being inhabited by a parallel spin, which 
occurs at Np = N s /2 sites. 

We introduce the possibility that the spin-blocked path can 
be opened by any process that alters the relative spin orienta- 
tion between the two sites. The probability for the blockade 
to be lifted by the time the next hopping attempt takes place, 
Xh, is p. We thus modify the effective density of P sites to be 
[1 — p]Np. The resulting modification of the density of sites 
that can be hopped to, N e /f, within the model of Miller and 
Abrahams, accounts for the possibility of spin flips of spins 
located at singly occupied sites. Using the effective site den- 
sity, we write the bonding criterion as 

B c = Ua 3 y 3 c (N - N P ) + 4%a 3 N P F y 2 pdy, (1) 
3 Jo 

where y c = r c /a is the dimensionless critical length which dic- 
tates the threshold resistance R c = Roe 2yc \ %h = VQ l e 2y is me 
hopping time. A quantity y CQ = (3B c /4%a 3 N) 1 / 3 is defined as 
the critical inter-site spacing in the absence of all spin effects. 
In general, y c cannot be isolated in Eq. ([T]) and the resultant 
MR can only be obtained numerically, however, in the dilute 
carrier regime (Np <C N), the MR obeys the analytic expres- 



sion 

l N P r y °i o 

MR ^2^-^ / y 2 [p(0)-p(H)]dy, (2) 

where y Cl = y Co (l — Np/N)~ 1 / 3 is the renormalized critical 
inter-site spacing. The MR scales linearly with the fraction of 
singly occupied sites. 

Most of the results reported below are based on the form 
of p that is appropriate if the spin flips in these organic ma- 
terials are caused by the hyperfine interaction (HI)(9j[l31; ex- 
pressions and implications appropriate for spin flips caused by 
the spin-orbit interaction are summarized briefly at the end of 
this Letter. Figure [2] emphasizes the main results of our the- 
ory. Panel (a) shows our calculations of MR for three differ- 
ent organic semiconductors: the small molecule TNF which 
has electrons as carriers (blue), the polymer derivative of PPV 
(orange) with hole carriers, and a generic material (black) that 
possesses a smaller localization length (1 A) than the other 
two (1 A< aTNF <appv as given in figure caption). This re- 
sult suggests that organic materials with small localization 
lengths yield the largest MR. Panel (b) contains our calcula- 
tions for different site concentrations and indicates that dilute 
site concentrations are more magnetoresistive. Panel (c) and 
(d) demonstrate the different MR structures and their depen- 
dence on the hyperfine field and site concentration obtainable 
from our theory. These results are derived and discussed in 
the remaining portion of this Letter. 

The case of fast hopping - When hopping is faster than the 
hyperfine precession frequency, the carrier spin experiences 
a random field for the short duration of time that it resides 
at a site. The spin-flip rate is identical to the well known spin 
relaxation rate from HI in the motional narrowing regime l25l : 



Xs V H+ X h 2Xh 

vh and v are precession frequencies due to the external field H 
and the in-plane internal hyperfine fields, of strength h. The 
probability for the P-spin to flip at the next hop is p(x s ) = 
1 — e~ Xh / Xs , which will permit a hop to the target site. This is 
a condition met in part when the density of sites is high. In 
this case the probability for a spin flip is p(i s ) w th/^s and the 
MR takes the form 

1 Np ry c i o 1 1 

MR^2^— / y 2 z h ( — — t-tW (4) 

y\ N Jo h \(0) x s (H) } y 

The hyperfine fields are random at each site so a correct de- 
scription of the MR involves an average over the Gaussian 
distribution of these fields. 

The resulting expression of Eq. ([4]) is cumbersome but can 
be simplified considerably by using the typical assumption 
that y Cx ^> 1 . The resulting MR response, averaged over the 
Gaussian distribution of hyperfine fields with width hi, is 

(MR) = \^vW c [1 - J- ln(l + v 2 H z 2 c)] , (5) 
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FIG. 2. (a) Calculated magnetoresistance % for a polymer PPV 
derivative (orange line) with localization length a = 3 A |24|, the 
small molecule TNF (blue line) with localization length a = 1.8 
A and a generic organic (black line) with localization length 
a = 1 A. Total site density is N = 10 26 m -3 and singly occupied 
site density is Np = 10 25 m -3 . (b) Calculated magnetoresistance for 
generic organic (black) and TNF (blue) at site densities N = 10 26 
m -3 (solid) and N = 2 x 10 26 m -3 (dashed), (c) Normalized magne- 
toresistance for generic organic at hyperfine field distribution widths 
(hj) of 0.5 mT (blue line), 1 mT (red line), and 2 mT (orange line). 
Site densities are same as in (a), (d) Normalized magnetoresistance 
for generic organic at total site densities of N = 5 x 10 28 m -3 (solid 
line) and N = 1 x 10 29 m -3 (dotted line). Only hyperfine interaction 
is assumed; unless noted otherwise, all figures use hj = 1 mT and 
v = 10 12 s" 1 . 



where l c is the hopping time at the critical radius and v/ is 
the precession frequency corresponding to the field hj. The 
positive MR can be understood by considering the field de- 
pendence of the relaxation mechanism; e.g. an increasing field 
suppresses spin relaxation via HI which makes the spin block- 
ade more effective. The dependence of MR on magnetic field 
here is identical to an earlier calculation for amorphous semi- 
conductors performed in the fast hopping limit l26l . Such 
agreement suggests our model provides an accurate descrip- 
tion of spin relaxation induced MR. 

The case of slow hopping - The condition 1 /v/ ^> Xh may 
not always be suitable for organic systems since the mobilities 
are so low. In positionally disordered systems the low mobil- 
ity is entailed by low site concentrations. As the site concen- 
tration is reduced, the site separations increase and the carrier 
hopping rate is reduced, leading to the condition 1/v/ < T c . 
During the requisite waiting time to hop, the carrier spin at i 
and target- site spin at j experience the applied field and their 
respective hyperfine fields hi and hj. Given two spins initially 
P aligned, the different hyperfine fields at the two sites ro- 
tate the spins to produce the possibility of AP alignment. We 
interpret this as a spin flip at either site. The time-averaged 
probability that the next hop is successful is p(H) = pij + pji 
where pij is the probability for the carrier at site i to be op- 
posite its initial state while the carrier at site j remains in its 
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The second term, pji is the reverse possibility. The MR, from 
Eq. ([2]), is ascertained to be 
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We note that the MR is independent of the hopping rate, which 
contrasts starkly with the fast hopping case. We now explore 
the MR lineshapes for these two cases. 

The magnetoresistance lineshape - Two characteristic fea- 
tures to quantify MR are its value at high magnetic fields, 
MR sat , and its width, 8, which we define as the half-width at 
half-maximum. The cross-over from slow to fast hopping can 
clearly be seen by viewing the saturated MR in Figure [3] and 
the MR width in the inset of the same figure. Slow hopping 
which results here from large inter-site distances, is conducive 
to large values of MR. The formula for MR sat in the low site 
density regime, derived from Eq. (|7|, is 
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Remarkably, the saturated MR is independent of HI. We note 
that in the deuterated PPV experiments of Ref. [29j 8 of 
the magnetoluminesence decreased when the hyperfine field 
was reduced whereas the high field magnetoluminesence was 
nearly unchanged. This result is consistent with our theory. 
In the fast hopping limit, the site concentration dependence 
is even stronger since MR sat contains an exponential depen- 
dence on N through x c . The MR widths also take on very 
different behaviors which are discerned from their MR ex- 
pressions above and shown in the inset of Figure [3] In the 
fast hopping region, the width is independent of the hyperfine 
interaction but strongly dependent on the hopping rate. Note 
that the width for fast hopping is larger - this is due to the 
quicker hopping rate which results in greater fields being re- 
quired to suppress HI spin relaxation. MR widths as large as 
40 mT - much greater than the hyperfine fields present - have 
been measured USED); Fig. [3] suggests these large widths are 
related to the hopping rate and not the hyperfine field. Often 
it has been assumed 1 29 , [31] [32) that 8 must depend on hi. We 
find this is only true in the slow hopping case where the MR 
curve roughly follows that of a Lorentzian of width hi. 

Spin-orbit coupling effects - For fast hopping, Eq. ^ is 
quite general in that any spin relaxation mechanism can be 
included for % s . Here we consider the influence of spin- 
orbit coupling (SOC); a recent calculation of SOC in sev- 
eral organic solids suggests that SOC may be significant l33l . 
Additionally SOC manifests itself by producing inhomoge- 
neous g-f actors (IG). The result of SOC produces a spin re- 
laxation rate: X7 1 = [^ 2 v^(v^ + x^ 2 ) _1 J r f f]\ l where y de- 
termines the SOC strength, b = ^3/108^, and 8g is known 
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FIG. 3. Percentage saturated magnetoresistance versus site concen- 
tration. Blue line: slow hoppingusing Eq. |8]); red line: fast hopping 
using high field limit of Eq. pi with hi = 1 mT, vo = 10 12 s _1 , 
a = 1 A, and Np = 10 25 m -3 . Inset: Magnetoresistance width, 5, as 
a function of site concentration for same parameters. 



to be proportional to y |25|. The field-independent portion 
reduces the total MR by but leads to no other qualitative 
change. This result is expected since SOC is field-independent 
and has been observed in Alq3 doped with Iridium [9 ]. IG, 
which increases with increasing field, leads to negative MR: 
-±^p 2 ln(l +v^x 2 ) where (v/ + v#)t c <C 1. At low fields 
this effect is expected to be small compared to Hi-induced 
MR. Recently, IG was studied in the slow hopping regime 
l34l but 8g was found to be unrealistically high to explain 
the magnetic field effects fT5ll . 

The prediction of large MR in the slow hopping regime ne- 
cessitates measurements over a controlled and wide range of 
site concentrations. We now discuss experimental strategies to 
observe our theoretical predictions. For r-percolation theory 
to be valid, the organic system must possess a low density of 
molecular sites and this density must be controllable. Conduc- 
tion via r-percolation was identified in TNF films by measur- 
ing the electron mobility through time-of-flight experiments 
lf2~Tl l22ll . The molecular density of TNF could be carefully 
controlled by dispersing TNF in a polyester host that did not 
alter the transport properties. We suggest similar experiments 
to look for the magnetoresistive dependences on hopping rate 
described in this Letter. We conclude by noting that this the- 
ory has implications for MR effects in amorphous semicon- 
ductors (25l and colloidal quantum dots |35], as well as for 
spin diffusion in organic spin valves 1 36 1 . 
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